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Abstract 

We consider parametrized families of linear retarded functional differential equa- 
tions (RFDEs) projected onto finite-dimensional invariant manifolds, and address the 
question of versality of the resulting parametrized family of linear ordinary differential 
equations. A sufficient criterion for versality is given in terms of readily computable 
quantities. In the case where the unfolding is not versal, we show how to construct 
a perturbation of the original linear RFDE (in terms of delay differential operators) 
whose finite-dimensional projection generates a versal unfolding. We illustrate the 
theory with several examples, and comment on the applicability of these results to 
bifurcation analyses of nonlinear RFDEs. 
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1 Introduction 



Differential equations are used to model a very wide variety of phenomena. Frequently, these 
differential models contain several parameters which are often varied or "tuned" to describe 
more accurately the phenomenon under study. Thus, there is considerable interest, from 
both a pure and an applied point of view, to understand how the properties of solutions of 
a parametrized family of differential equations are affected by variation of the parameters. 
This philosophy is at the core of bifurcation theory. 

One large class of differential equations which are particularly important in applications 
are retarded functional differential equations (RFDEs) [12|, [13[ , which includes the class of 
ordinary differential equations (ODEs), the class of delay differential equations, as well as 
certain types of integro-differential equations, among others. These equations are used to 
model various phenomena in fields ranging from mathematical biology |Q, [II], [T5[ [T(| to 
industrial processes [fT8 1 , and to atmospheric science Rl9| . The theory for both linear and 



nonlinear RFDEs is rather well-developped []12|, |13|j . Essentially these equations behave like 
abstract ODEs on an infinite-dimensional (Banach) phase space. Thus, many results which 
are known for ODEs on finite-dimensional spaces have analogs in the context of RFDEs. For 
example, in the neighborhood of an equilibrium point of a nonlinear RFDE, there exists local 
invariant manifolds (stable, unstable and center manifolds) which are tangent to the corre- 
sponding invariant subspaces of the linearized equations about the equilibrium point, and 
on which the flow near the equilibrium is either exponentially attracting (stable manifold), 
exponentially repelling (unstable manifold), or non-hyperbolic (center manifold). In the con- 
text of bifurcation theory, the center manifold reduction of the flow is important, since this 
is where bifurcations (qualitative changes in the flow) take place as parameters are varied. 
In applications, there have been many studies (see for example 0) of specific RFDEs where 
stability and bifurcations of equilibria were investigated using the center manifold reduction 
theory developped in |I2j and [Tj|. Another example of the similarity between ODEs and 
RFDEs is in || [7|, where the theory of Poincare-Birkhoff normal forms was extended to 
RFDEs. 

One aspect of the bifurcation theory of RFDEs which, surprisingly, has not yet been 
developped is that of extending Arnold's theory |IJ on versal unfoldings of matrices to the 
case of parameter-dependent linear RFDEs. This is the purpose of our present paper. Versal 
unfoldings of RFDEs for certain singularities have been computed for particular classes of 
RFDEs in the study of restrictions on the possible flows on a center manifold see [||, Bl [T7| . 
These unfoldings are computed using the normal form theory of Faria and Magalhaes JF]. 
However no attempt is made to give a systematic treatment of unfolding of linear RFDEs. 

We adopt the following strategy. We begin with a linear RFDE x(t) = £o(x t ) (where 
£q is a bounded linear functional operator) whose semiflow restricted to a finite-dimensional 
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subspace is defined by the matrix B. Using a versal unfolding of B, we explicitly construct 
a parametrized family C(a) of bounded linear functional operators whose finite-dimensional 
restricted semiflow is defined by a versal unfolding 13(a) of the matrix B. In comparison, the 
realisation of linear ODEs by linear RFDEs obtained by Faria and Magalhaes || provides 
an existence result. Using the Hahn-Banach theorem they show the following. For any finite 
dimensional matrix B, a necessary and sufficient condition for the existence of a bounded 
linear operator £ from C([— r, 0], IR n ) into M. n with infinitesimal generator having spectrum 
containing the spectrum of B is that n be larger than or equal to the largest number of 
Jordan blocks associated with each eigenvalue of B. While some of the proofs in our paper 
have a flavor similar to the ones in ||, our results cannot be deduced from realisability 
results. 

Recall that for a given cxc matrix B with complex entries, a p-parameter unfolding of 
B is a p-parameter analytic family of matrices B(a) such that B(ao) = B for some cto G C p . 
The unfolding B(a) is said to be a versal unfolding of B if for all g-parameter unfoldings 
A(/3) of B (with A(/3q) = B), there exists an analytic mapping : C q — > C p and an analytic 
family of invertible matrices C(f3) satisfying 

0(A)) = «o 

C(A0 = / 
A((3)=C((3)B(m)(C((3))- 1 . 

Thus, a versal unfolding of B is, up to similarity transformations, a general analytic pertur- 
bation of B. A versal unfolding is said to be mini-versal if the dimension of the parameter 
space is the smallest possible for a versal unfolding. Of course, the concept of versal unfold- 
ing of matrices is of importance not only to linear differential equations, but to nonlinear 
differential equations as well, since questions of stability and bifurcations of equilibria in 
nonlinear systems always involve an analysis of an associated linearized system. It is there- 
fore important to understand the dependence of the associated matrix on system parameters 
(i.e. in the case where the resulting unfolding is not versal, there are restrictions on the 
movement of the eigenvalues, which may influence the possible range of dynamics). 

In the space Mat cxc of c x c matrices with complex entries, let £ denote the similarity 
orbit of the matrix B. We will use the following sufficient criterion for versality which can 
be found in and @]. 

Proposition 1.1 Let B(a) be a p-parameter unfolding of the matrix B. If 

Mat cxc = TS B(ao) + D a B(a ) ■ C p , 
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where «o is such that B(olq) = B and TXg( Q0 ) denotes the tangent space to S at B(a>o), 
then B(a) is a versal unfolding of B. If the codimension of THB( ao ) i s equal to p, then the 
unfolding is mini-versal. 

As motivation for our analysis in this paper, consider a parameter dependent nonlinear 
RFDE 

x(t) = F(x t ,a) (1.1) 

which has a non-hyperbolic equilibrium, which we assume for the sake of simplicity to be at 
x = 0, a = 0. As previously mentioned, a crucial step in the analysis of the bifurcations and 
stability of this equilibrium is to perform a parameter dependent center manifold reduction 
of the RFDE J7| in order to obtain a parameter dependent nonlinear c-dimensional ODE 

z(t) = Bz + G(z,a), (1.2) 

where B is a cx c constant matrix, and G is nonlinear in z and a. The emphasis is important, 
since terms which are parameter dependent and linear in z are contained in the expression 
for G and they are the ones which unfold the matrix B. Thus, it is natural to investigate 
the versality of this unfolding, for reasons that we have previously mentioned. However, 
the potential difficulty here is that the reduction process for RFDEs to finite-dimensional 
invariant manifolds introduces restrictions on the possible types of nonlinearities G which 
can be achieved in ( p..2|) || . Thus, there is no a priori guarantee that a versal unfolding of 
the matrix B can be achieved in G by this reduction process, even if we have perhaps many 
parameters (e.g. more than the codimension of the singular matrix) in the original RFDE. 
It is clear that only terms which are linear in x t in the right-hand side of ( |1 . 1| ) contribute to 
terms which are linear in z in the right-hand side of (|1.2|) . Thus, we restrict our attention 
to the case where ( |1.1|) is a parametrized family of linear RFDEs. 

We have two main results. The first gives a sufficient condition on the parametrized 
linear RFDE (|1.1| ) which guarantees that the right-hand side of (|1.2j ) is a versal unfolding of 



the matrix B. The second main result is twofold: first, we show that despite the previously 
mentioned restrictions on G in ( p..2|) , it is always possible to realize a versal unfolding of 
B by a suitable choice of parameter-dependent RFDE in ( |1 . 1| ) ; and furthermore, we give 
a "canonical" method of computing such a parameter-dependent RFDE in terms of delay 
differential operators. The theory is then illustrated with several examples. 

Although our original interest and motivation lies in bifurcation theory (i.e. B in ( |1.2j ) 
has all of its eigenvalues on the imaginary axis), there is no additional complication in 
considering a reduction of ( |1 . 1| ) to a general finite-dimensional invariant manifold (i.e. not 
necessarily a center manifold). Therefore, we develop the theory in this general context. 

As is the case for Arnold's theory of versal unfoldings of matrices, we develop our theory 
by working in complex spaces; since the diagonalization theory is much simpler in this 
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context. Versal unfoldings in the real case can be constructed by a decomplexification of the 
complex unfolding, as is done in |IJ and which we illustrate in Section 7. 



2 Reduction of linear RFDEs 

Let C n = C([—t, 0], C n ) be the Banach space of continuous functions from the interval [— r, 0], 
into C n (r > 0) endowed with uniform norm. We are interested in the linear homogeneous 
RFDE 

z(t) = Co(zt), (2.1) 
where Co is a bounded linear operator from C n into C n . We write 

Coifp) = f dri(e)<p(o) t 



where rj is an n x n matrix- valued function of bounded variation defined on [— r, 0] . Let 
A denote the infinitesimal generator of the semiflow generated by equation ( |2.1|) . Then 
it is well-known that the spectrum cr(A ) of A is equal to the point spectrum of A , and 
A G cr(A ) if and only if A satisfies the characteristic equation 

detA(A)=0, where A(A) = XI n - / dr](9)e xe , (2.2) 

where I n is the n x n identity matrix. We suppose that A C C is a non-empty finite set of 
eigenvalues of A , with corresponding generalized c-dimensional eigenspace P. Using adjoint 
theory, it is known that we can write 

C n = P@Q (2.3) 



where Q is invariant under the semiflow of (2.1), and invariant under A$. 



Define C* = C([0, r], C n *), where C n * is the n-dimensional space of row vectors. We have 
the adjoint bilinear form on C* x C n : 

v)n = mm - y° J ^(e - 0)<m*mo#- ( 2 - 4 ) 

We let $ = (ipi, . . . , <f c ) be a basis for P, and ^ = col^i, • • • , V'c) De a basis for the dual 
space P* in C*, chosen so that (\1/, $) n is the c x c identity matrix, I c . In this case, we have 
Q = {tp G C n : (^j V 9 )™ — 0}. We denote by i? the c x c constant matrix such that $ = $£>. 
The spectrum of B coincides with A. Using the decomposition ( |2.3j ), any z G C n can be 
written as z = $ x + y, where igC and y G Q is a C 1 function. The dynamics of (|2.1| ) on 
P are then given by 

x = Bx. 
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3 Parametrized families of linear RFDEs 



Consider now a smoothly parametrized family of linear RFDEs of the form 

z{t)=C{a){zt), (3.1) 

where a G C p , and C(ao) = Co is as in (j2.1|), for some a G C p . In the sequel, we will assume 
that a translation has been performed in the parameter space C p such that «o = 0. Our 
outline and notation here follows closely that of [0. We rewrite as the system 

z(t) = Co(zt) + [C(a)-Co](zt) 

(3.2) 

d(t) = 0. 

The solutions of this system are of the form z(t) = (z(t), a{t)) T G C n+P (where the superscript 
T denotes transpose), the phase space is C = C n+P = C([—r, 0], C n+P ), and we write ( |3.2|) as 

k{t) = C z t + F{z t ), (3.3) 

where C ((u,v) T ) = (C (u),0) T and F((u,v) T ) = ([£(v(0)) - C )(u), 0) T , u G C n , v G C p . 

Let A, P, Q, $, ( , ) n and 5 be as in the previous section. Define P = P x C p , 
Q = Q x R, where R = {v G C p : f(0) = 0}, and consider for bases of P and P*, 
respectively, the columns of the matrix $ and the rows of the matrix 

• =i * ), *=f* 

I v ) \ 7 P 

which satisfy (\1/, $) n +p = Ic+p- We have $ = $S, where I? = diag(5, 0). It follows that we 
have an invariant splitting C = P © Q. 

Let BC n denote the space of functions from [— r, 0] to C" which are uniformly continuous 
on [— r, 0) and with a jump discontinuity at 0. If we define X : [— r, 0] — > C n by 



X o {0) 



In = 

-r<9<0, 



then the elements of BC n can be written as £ = if + Xo/i, with </? G C n = C([— r, 0],C") 
and /i G C™, so that -BC n is identified with C n x C™. In order to study ( |3.3|) , we need 
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to consider the space BC = BC n x BC P , which can be identified with C x C n+P . Define 
Y : [-r,0] -^C?by 



Y o (0) 



I p 6 = 



o -t <e <o, 



Let 7r : 5C n — > P denote the projection 

7r(^ + Xofl) = <^) n + *(0)/i], 

where G C n and /i G C". We consider the projection 7r : BC — > P given by 



*(0) 



n+p 



/' 

V 



Tl((p + X fJL) 

^(0) + z/ 



We now decompose z in ( |3.3|) according to the splitting 

BC = P® ker n, 
with the property that Q C ker7f, and get 



.r 
a 

V 
w 



B 



X 


+ ^(0)F ^<f> 
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a 






w 





w 



+ (/-7f) 



*0 



.r 

a: 



+ 
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(3.4) 



where x G C c , a G C p , y G Q 1 = Q fl C£, w G R 1 = R f] C p (C£ and C£ denote respectively 
the subsets of C n and C p consisting of continuously differentiable functions), and Aqi is the 
operator from Q 1 = Q (1 C 1 = Q 1 x R 1 into ker tt defined by 



Ar>i 



+ 



X 
Y 



Co 



0(0) 
V>(0) 



If Aqi : Q 1 C ker 7r — ► ker7r is defined by Aqup = <fi + X [£ if) — 0(0)]; then ( |3.4| ) is 
equivalent to 











d 




d 


1/ 




eft 







Bx 




+ 



*(0)[£(a(0) + w(0)) - Co]($ x + y) 




A-Q^y 
w - Y w(0) 



+ 



(/ - tt)X [£(a(0) + iu(0)) - C ]($x + y) 
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Since w G R, it follows that w(0) = 0, so that we get the following equations in BC n = 
P © ker 7r 

i = Bx + ^(0)[£(a) -£ ){$x + y) 

d (3 ' 5) 
—y = A Q xy+(I-ir)X [£(a)-£ )(®x + y), 

where x G C c and y E Q 1 . 

4 Reduction to parameter dependent invariant mani- 
fold 



In this section, we show that fl3.4| ) admits a local, semiflow invariant, c + p-dimensional 
manifold in BC, which is tangent at the origin to P, and such that the dynamics of ( |3.4| ) 
restricted to this manifold are linear in x G C c . 

We want the nontrivial part of our invariant manifold to be of the form 

u = $x + h(a)x (4.1) 

where h : C p — ► Mati xc (<5 1 ) is a smooth map, and MatixclQ 1 ) denotes the space of 1 x c 
matrices whose elements are in the space Q 1 which has been defined in the previous section. 
Now, as an infinite-dimensional system, the RFDE can be written as || 

on 

— =Au + X [£(a)-£ ]u. (4.2) 



Combining (Op and the first equation in Q3.5| ) we obtain 

"i = + (4.3) 
= ($ + h(a))(Bx + V(0)[£(a) - £ ]($ + h(a))x). 

The expression (|4.1|) represents a locally semiflow invariant manifold of (|3.2|) near the origin 



z, a) = (0, 0) if equation ( |4.2|) is equal to ( |4.3|) for values of a in a small neighborhood of 0. 



Rearranging the two equations and simplifying implies that we need to solve 
(A Q1 7t)X [£(a) - £ ])h(a) + (I - n)X [£(a) - £„]($) 

(4.4) 

-h(a)(B + #(0)[£(a) - £ Q ]($ + /i(a))) = 
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for small values of a in a neighborhood of 0. 

Our main tool for proving this result is the implicit function theorem (IFT). Because of 
the smoothness properties required for the application of the IFT, we need to work with the 
C 1 norm instead of the uniform norm on the space Q l . We claim that the spectral properties 
of the operators Aq, A and Aqx which are given in lemmas 5.1 and 5.2 of || remain valid if 
we replace the uniform norm by the C 1 norm in the domain of these operators. In particular, 
we have 

a(A Q ,)=Pa(A Q ,)=a(A )\A, (4.5) 

where a denotes the spectrum, and Pa is the point spectrum. Clearly the point spectrum 
does not depend on the norm. Also, the resolvent set does not change, since the above 
operators are bounded when we replace the uniform norm by the C 1 norm in the domain of 
the operators, and for every A in the resolvent sets (with respect to the uniform norm) of 
these operators, the resolvent operator (A — A/) -1 (where A is any of these operators) is 
defined on all of BC n , and not just a dense proper subset. Therefore, with respect to the C 1 
norm, A — A / is bounded and surjective, so (^4 — A J) -1 is bounded, by Banach's theorem. 

Proposition 4.1 Let (Q l ) c = Q 1 x ■ • • x Q 1 (c times) endowed with C 1 norm 

c 

\\(h\...,h c )\\ mc = Y,(\h% + Mc), 

i=l 

(where \ \c denotes uniform norm), and let (ker7r) c = ker7r x ■ ■ • x ker7r (c times) endowed 
with norm 

c 

IK^ 1 + X a\ . . . ,r + X a c )|| (ker7r)c = (l^lc + l«V). 

8=1 

Consider the nonlinear operator N : C p x (Q 1 ) — > (ker7r) c defined by 

N(a,h) = (A Q i + (I-Tr)Xo[C(a)-£ ])h + (I--K)X [C(a)-C ]($)- 

(4.6) 

h (B + *(0) [£(<*) - Cq]{$ + h)) , 

where the actions of Aqi and C(a) on h e {Q l ) c are defined componentwise in the obvious 
way. Then, there exists a neighborhood V of in C p and a unique smooth mapping a i — ► 
h(a) from V into (Q 1 ) such that h(0) = and such that N(a, h(a)) = for all a 6 V. 

Proof It is clear that with the chosen topologies, iV is a smooth mapping. Moreover, it 
is also clear that iV(0,0) = 0. The partial Frechet derivative ^(0,0) is the bounded linear 
operator from {Q 1 ) into (ker7r) c defined by 

N h (0, 0)v = Jv = A Q iv - vB. 
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We suppose that we have chosen a system of coordinates such that the c x c matrix B is in 
Jordan canonical form 

/ A x a, \ 

A 2 <7 2 



B 



V 

where a,i = 1 or cr, = 0. Suppose that v 
Then this is equivalent to 



•• cr c -i 

v c ) G (Q 1 ) c \ {0} is such that Jv = 0. 



A Q i{v x , . . . ,v c ) = (X 1 v 1 ,a 1 v 1 + A 2 w 2 , . . • , a c _i^ c 1 + X c v c ), 

which implies that one of the Aj must be in the point spectrum of Aqi , which is a contradiction 
(see ([Of)). Therefore, ker J = {0}. 

We now show that J is surjective. We use an approach similar to the proof of Theorem 
5.4 of p. Let i G {l,...,c}, <p G Q 1 , and £ G ker7r be such that £ = (Aqi — Xil)ip. 
Define v = (v 1 , . . . ,v c ) G (Q 1 ) by v k = 0, k ^ i, and t>* = ip. Then f = Jv has the form 
/ = (/V--,/ c ), where 



f if A; = i 
f k = { —<Ji(f if k = i + 1 and i < c 
otherwise. 



(4.7) 



Let g = (g 1 , . . . , g c ) G (ker 7r) c . We now use an induction argument to show that there exists 
an h = (/t 1 , . . . , h c ) G (Q 1 ) which is such that Jh = g. We know from Q4.5p that Ai is in the 
resolvent set p(Agi). Therefore, there exists ip\ G Q 1 such that [Aqi — Xil)<pi = g 1 . If we 
define H x = (<pi,0, ■ ■ ■ ,0) G (Q 1 )", then we get from (fT7]) that (Jfti) 1 = g 1 . Suppose now 
that i G {1, . . . , c — 1} is such that there exists G (Q 1 ) c satisfying 



; JHA P = g p , \/p < i. 



Define H 4+1 = (JTtiY +1 e ker7r. Since X i+ i G p(Agi), there exists <£> i+1 G Q 1 such that 
(A Q i - X i+1 i> i+1 = g i+1 - H i+1 . If we define G l+1 = ((\ . . . , C c ) G (Q 1 ^ by = if 
k^i + 1, C +1 



it follows from ( |4.7| ) that 

gi+l_ H i + l ifk=i+l 



(JG, 



i+lj 



-(Jitpi + i if k — i + 2 and z + 1 < c 
otherwise. 
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If we now set TLi+i = Hi + Gi+i, we get from 

gL+i Hp — %-\-\ 



g p if p < i. 



By induction on i, we have that (fOf) holds for i = c, and thus J is a surjection. It now 
follows that J has a bounded inverse, and we get the conclusion of the Proposition by virtue 
of the implicit function thereom. I 

Proposition 4.2 Let h(a) be the solution of N = in \j.b\ j defined for all a G V , where V 
is some neighborhood of the origin in C p . Then the following set 

W = {((x, a), (y, w)) G C c+P x Q 1 : w = 0, y = h(a)x, x G C c , a G V} (4.9) 

is a locally semiflow invariant manifold for the system ( \3.4\ )- The (non-trivial) dynamics on 
this manifold reduce to 

± = Bx + (0)[£(a) - £ ]($ + h(a))x. (4.10) 

Proof It is obvious that W is tangent to P. The semiflow invariance of W follows from 
substitution of y = h(a)x into ( p.5| ), and using the fact that h(a) is a solution to N = in 
L6|). ■ 



5 Sufficient condition for versality 

From the previous section, we have seen that the dynamics of ( |3.1| ) near the equilibrium 
solution (z, a) = (0, 0) reduces to the c-dimensional parametrized linear system 

x = B(a)x, (5.1) 

where by ( |4.1C|) we have 

B(a) = B + *(0)[£(a) - £ Q ]($ + h(a)). (5.2) 
Note that B(0) = B and since h(0) = 0, we have 

D a B(0) = 9(0) D a [£(a)(<l>)}\ a=0 . 

Let Mat cxc denote the space of c x c matrices with complex entries. For i, j G {1, . . . , c}, 
let G Mat cxc be the matrix whose elements are all except the element in row i and 
column j, whose value is 1. For Ui, U 2 G Mat cxc , denote [Ui, U 2 ] = UiU 2 — U 2 U\. 

Denote the linear mapping 6 : Mat cxc — > C c by Q(Eij) = e^-i) c+ j, where ee denotes 
the row vector whose components are all except the £ th which is 1. 
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Theorem 5.1 The parametrized family \3. J\ ) generates a versal unfolding B(a) (see (\5.^{ )) 
of the matrix B = (\I>, $) n if the (c 2 + p) x c 2 matrix 



( 



S 



\ 



e([B,E 12 ]) 
&([B,E CC }) 

e( *(o)^- [£{<*){*)] 



e(*(o) A [£(«)($)] 



a=0 



a=0 



/ 



/ias ranA; c 2 . If, in addition, 

dim span(9([ J B, E u ]), Q{[B, E l2 }), ...,G([B, E cc ])) = c 2 - p, 
then the versal unfolding 13(a) is mini-versal. 

Proof In Mat cxc , the tangent space of the similarity orbit through B is given by |L] 

= {[B,X] : XeMaU}. 
Thus, if the matrix S has rank c 2 , we have 

Mat cxc = O + (0) D a [£(a)($)]\ a=0 ■ C p , 



(5.3) 



which implies that the mapping a i — > £>(a) is transversal (mini-transversal if ( |5.3|) holds) 
to the similarity orbit of B at a = 0. The conclusion now follows from Proposition |1.1| . I 

From the previous result, we are now led to define a notion of versal unfolding for RFDEs 
of the type JO). 



Definition 5.2 The parametrized family of RFDEs \3. 1\ ) is said to be a A-versal unfold- 
ing (respectively a A-mini- versal unfolding,) for the RFDE ($.J\) if the matrix B(at) defined 
by (\5.3J is a versal unfolding (respectively a mini-versal unfolding) for B. 
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6 Decomposition of Mat cXc by \I/(0) 

Consider now the following problem: given a linear homogeneous RFDE such as ( |2.1|) and a 
set A of solutions to the characteristic equation (|2.2|) , find a A-versal unfolding flOl ) for (|2.1|) . 
It is certainly not immediately obvious that this problem need admit a solution, since the 
structure of the right-hand side of (|4.10|) is severely restricted by the structure of the matrix 



^(0). To solve this problem, we will need to characterize the subspace of matrices in Mat cxc 
whose columns are in the range of ^(0), and show that one can build a versal unfolding of 
the matrix B in ( |4.10[ ) even in this restricted context. We start with the following 

Definition 6.1 We define ^(^(O)) to be the set of all c x c matrices whose columns are in 
the range of the matrix ^(O). 

The main result we will need in order to construct the above-mentioned versal unfolding is 
the following: 

Proposition 6.2 Let T : Mat cxc — > Mat cxc be defined by T(M) = [B, M) = BM - MB. 

There exists a subspace W C 7^.(^(0)) such that 

Mat cxc = range(T) © VV. (6.1) 



This entire section is devoted to proving this result, including giving an explicit construction 
of W. We will start by first proving several lemmas, and then proceed to the proof of 
Proposition |6.2| . 



6.1 Commutator of B 

We assume that the matrix B is in the following Jordan block diagonal form 

B = diag(5 1 1 (A 1 ), • • • , <(A!), S?(A a ), . . . , 5 fc 2 2 (A 2 ), . . . , B[(X r ), . . . , B r kr (X r )), (6.2) 

where Ai, . . . , A r are the distinct (without multiplicities) eigenvalues of B, and the Jordan 
block Bj(Xj) is n jti x n jti of the form \j I njtt + N nuXn . e , where N n . lXnj l is the nilpotent 
matrix with l's on the upper diagonal, and 0's everywhere else. Moreover, we assume that 
for each j 6 {1, . . . , r}, we have 

%i > n jt2 > ■ ■ ■ > n j)kj . 



The first result we need is the following lemma, which can be found in [10] 
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Lemma 6.3 Let B* denote the conjugate transpose of B, and let M e Mat cxc be such that 
[B*, M] = 0. Then M is of the form 

M = diag(.M At), (6.3) 

where Aij is (n^i + • • • + rij^.) x (n^i + • • • + Uj^ ) matrix, partitioned into blocks with 
dimensions Tlj p /n Tlj and of the form illustrated in Figure [Z], where each oblique segment 
in each separate block denotes a sequence of equal entries, and all other entries are zero. 

n. , n. „ n. „ n. , 

J,l J,2 j,3 j,4 



V 
V 

V 

Figure 1: Structure of a matrix Aij. In this example, kj = 4, n^i = 9, rijp = 5, nj,3 = 4, 
rij-,4 = 2. 




The second lemma we need is the following. The proof is left to the appendix. 



Lemma 6.4 Let B be as in (\6. 2 ), and consider the mapping T : Mat c> 
by T(M) = [B, M] . Then Y 6 range(T) if and only if Y is of the form 



Y 



( yi,2 ■■■ 

3^2,1 3^2,2 " " " 

\ y r ,i y r ,2 ■ ■ ■ 
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3V 



\ 



■ ■ ■ y r , r J 



Mat cxc given 



(6.4) 



where y pA is a (n Pt x + • • • + n p ± ) x (n 9i i + • • • + n q ^ ) matrix, with the only constraint being 
that for each j e {1, . . . , r}, the sub-matrix is partitioned exactly as Aij (see Figure 
and is such that the sum of the elements in each given oblique segment is zero (however, in 
contrast to Aij, the elements o/D^i,j which are not on the oblique segments are completely 
arbitrary, i.e. not necessarily zero). 

We now define useful integers which give the number of columns of the matrix B correspond- 
ing to the first j eigenvalues. Let N = 0, and 

kj 



Using Lemma |6.4j , we now define a set of matrices S which forms a basis for range (T). The 
set S consists of all c x c matrices Y which are partitioned as in ( |6.4j ) and whose elements 
are all zero except one element whose value is 1, and possibly one other element whose value 
is —1, satisfying the following constraints: 

• if the element whose value is 1 is in the block where k ^ £, then it is the only 
non-zero element in the matrix Y, 

• if the element whose value is 1 lies in the block 3^,j for some j G {1, . . . , r}, then: 

— if this element (whose value is 1) is not on any of the oblique segments of Figure 
[I], then it is the only non-zero element in the matrix Y 

— if this element (whose value is 1) does lie on one of the oblique segments of Figure 
[I], then it does not lie in any of the following rows 

Nj-t + rij,!, Xj i • • ii h2 - Nj-x + n^i -I \-n j)kj , (6.5) 

and there is a ' — 1' at the bottom end of the oblique segment which contains the 
T. 

Note that the row numbers given by ( |6.5|) correspond to the last row of rij^. x nj^. blocks 
as shown in Figure [I] for each j . 

Now consider W C Mat cxc , where W is the subspace of matrices of the form 

diag(u;i, . . .,w r ), 

where Uj is a (n.,-,1 + • ■ ■ + %,&■) x (rij t x + • • • + ) matrix, partitioned into blocks with 
dimensions n^ v x rij jq , and of the form illustrated in Figure |^. In this figure, the only elements 
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which are not forced to be zero are those at the bottom of each oblique segment, illustrated 
with a bold dot. It is well-known [|1(| that 



r kj 

dim(W) = Yl ( 2t ~ 1 ) n ^" ( 6 ' 6 ) 

3=1 1=1 

n. , n. „ n. „ n. , 

J.1 J.2 j,3 j,4 



V 



n j,3 
V 

Figure 2: Structure of a matrix ujj. In this example, kj = 4. The only elements which are 
not forced to be zero are those illustrated with a bold dot, at the bottom of each oblique 
segment. 

It follows from the definition of the basis S of range(T) that 
Lemma 6.5 

Mat cxc = range(T) © W. (6.7) 
6.2 Characterizing the range of ^(0) 

Before we can prove Proposition we need to establish some properties about the subspace 
n(V(0)) defined in Definition |IJ 



























•> 


• • 
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Lemma 6.6 The matrix ^(O) = col(^i(0), . . . ,^ c (0)) satisfies the following property: for 
each j G {1, . . . , r}, if % G { + n j:1 , Nj-i + n j:1 + n j:2 , ■ ■ ■ , iV^-i + n^i H h ra^ } t/ien 

^(0)^0 and ^(0)A(A i ) = 0, 

where the characteristic matrix A(A) zs as in ( $.2J . Moreover, for each fixed j G {1, . . . , r} ; 
the set of row-vectors 

{ fpNj^+nj,! ^-1+^1+^,2(0), . . . , ip Nd _ 1+njA+ ... +njtk , (0) } 

is linearly independent in C n *. 

Proof Each row ipi{6) of the matrix ^(6) = e~ Be ^(0) must satisfy 

^(0) = - f ^(-9)^(9). (6.8) 



Since B has the form (O), it follows that for each j G {1, . . . , r}, if % G { iV 3 -_i ^j-i + 

n i,i + n j,2> • • • 5 ^j-i + n j.i + ' ' ' + n j,k } then the i th row of the matrix e~ is zero except for 
the diagonal element which is equal to e~ Ajfi '. It follows that for alH G { Nj-i + n^x, Nj_x + 
n i,i + n i,2> • • • ) Nj-i + n jtl H h n j:k . } we have 

i/>i(e) = e- x *%(0), 



which when substituted into (|6.8| ) yields ^(0)A(Aj) = 0. Recall that the columns of ^(0) 
form a basis of P*, therefore the set 

{ ^Nj^+n^^^Nj-t+n^+n^), ■■■> ^N j - 1 +n j , 1 +-+n jtk . {9) } 

is linearly independent in C([0, r], C n *), it follows that 

+Hj l+Jlj 2 

(0),...,Vat,_i +"j,iH l-^j.fc,- 

(0)} 

is linearly independent in C n *. ■ 
Remark 6.7 For each j G {1, . . . , r} ; /et 

{^JVj_i+^,i(0),^_i +fij i+rij 2 (0),...,^iv,_i +«j,lH hrij^fc . (0)} 

6e as m Lemma \6. b\ Consider the linear mapping 

Uj : C n — ► 

defined by 

(0),...,ViV f _i +n.j i iH hi^fc . 

(0)) ■ u. 

Then it follows from Lemma \6.b] that ILj is onto C kj . 
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6.3 Proof of Proposition p.'4 



We are now ready to prove Proposition |6.2| . For purposes of clarity, we first prove it in the 
case where r = 1 in ( |6.2|) , and then we show how to generalize the arguments of this proof 
to the case r > 1. Before giving the proof, we first establish some useful notation. 
In the case where r = 1 in (|6.2|), any matrix M G Mat cxc can be partitioned as 



/ M 1A Mi, 
M 2 ,i 



M 



2 

M 2 ,2 



M 2M 



\ Mk lt i M kl , 



Mk u k! 



(6.9) 



/ 



where M^,\ is n\£ x ni t x, for £, A G {1, . . . , k{\. It is convenient to label the elements of 
M according to this partitioning as follows: for £, A G {1, . . . , ki}, m G {1, . . . ,^i,a} and 
u G {1, . . . , we denote by M^' x,u,m the element of the matrix M which lies in the block 
M^x in (|6.9p , at the intersection of row u and column ni t \ — m + 1 relative to this block. 

Fix £, A G {1, . . . , ki} and let A) be the set of integers m G {1, . . . , ni } \} such that 
column ni y \ — m + 1 intersects an oblique line in the block M^\. In particular, 




ifni^ > n 



■ ni£ + 1, . . . ,rii jA } ifni )? < ni,A- 



Now, for all £, A G {1, . . . , fci} and m G Q(£, A), we define ft^ } x,m as the c x c matrix whose 
entries are all zero except ^\'x^' m = 1- It follows that the set of all matrices fl^x,m thus 
defined is a basis for the subspace W in ( |6.7| ) (see Figure 0). 

Similarly, fix £, A G {1, . . . , ki} and m G {1, . . . , Let V(£,\, m) consist of all 

u G {1, • • • , £ — 1} such that row u has an intersection with an oblique line in the block 
(£, A) at column ni } x — m + 1. Of course, if nj^ = 1, m G" <2(£, A) or m = min Q(£, A) then 
P(£,A,m)=0. 

For all £, A G {1, . . . ,ki}, choose m G <2(£, A) such that V(£, \,m) is nonempty and 
choose u G V(£, A,m). Define the c x c matrix E^,x,u,m with only two nonzero entries as 
follows. Let Ef'*'™'™ = 1, and ^'a,'™ 1 ^'^" 1 = where p(m,u) = m — ni£ + u. Note that 
p(m, u) is always less than m. See Figure |0| for an illustration of Q(£, A), "P(£, A, m) and 
p(m,u). Each of these matrices E^ x,u,m defined above is in range(T) (Lemma 



Remark 6.8 It follows from the definitions of the matrices ^ t x,m and the matrices E^,x,u,m 
that if A G {1, . . . , ki}, £ G {1, . . . , hi}, m G {1, . . . , n^x} and u G V(^, A, m) ; £/ien p(m, u) G 
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m = 6 




p(6,3) 



Q(£,A) 

Figure 3: Example of sets Q, P and of a point p(m, w) for a 5 x 7 matrix. Here Q 
{3,4,5,6,7} and V = {2,3,4}. 

Q(£, A) and i%,A,u,m + ^f,A,p(m,«) a matrix whose only non-zero element is 

( zp _i_ Q \£,A,u,m -i 



We now proceed to the proof of Proposition |6.2| in the case where r = 1. 



Lemma 6.9 Suppose B in ( \6.£\ ) is such that r — 1. Then Proposition \6.3j holds. 

Proof We begin by defining a subspace of matrices contained in 72.(^(0)) and isomorphic 
to W. We do this in the following way. 

By virtue of Remark 3/7, let t>i, . . . , G C n be such that Il 1 (^) is a fcx-dimensional 
vector whose only non-zero component is the £ th component, whose value is 1. Now, consider 
the c-dimensional vector \I/(0)t^. We label the components of this vector using the integers 
£ G {1, . . . , ki} and u G {1, . . . , ni^} as follows: (\I/(0)f£)^' u is the ■u*' 1 component of ^(0)^ 
if £ = 1; and the (n x i + • • • + Jii,£-i + w) <?i component of ty(0)v£ if £ > 1. We then have 



7c, u,e if u ^ {ni,i, . . . , rii,*^} 
(*(0) = < 1 if £ = £ and it = n u 



(6.10) 



otherwise 
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where the exact values of the coefficients 7£ lWr g are not important. 

Now for all £, A G {1, . . . , ki} and m G Q(£, A) define R^x, m to be the n x c matrix whose 
("-1,1 + • ■ ■ + ni,A — m + column is and all other columns are zero. We then define a 
linear mapping £ : W — > 72.(^(0)) by the following action on the basis elements of W: 

£(fi € ,A,m) = *(0)i2g,A,m, e,Ae{l,...,M, me Q(£,A). (6.11) 

It is clear that S is an isomorphism between W and W = £ (W) C 72.(^(0)), since the set 

: e,Ae{l,...,fci}, mGQ(£,A)} 

is linearly independent in Mat cxc . 
Our strategy now is to show that 

W C range (T) + VV, 

from which it follows from ( |6.7|) that 

Mat cxc — range(T) + W, 
and since W and W are isomorphic, we will then have 

Mat cxc = range(T) © VV. (6.12) 

Any matrix Z G Mat cxc can be written as a sum of two matrices as is illustrated in 
Figure ^, where the elements which are not on the oblique segments are all zero. It is clear 
from Lemma jA that the second summand in Figure [| belongs to range(T). Thus, for any 
Z G Mate xc, we define F(Z) to be the first summand in this decomposition, as illustrated in 
Figure |j. It immediately follows that for any Z G Mat CXC ! we have Z — T(Z) G range(T). 

Fix a value of A G {1, ... , fci}. For all £ G {1, . . . , fci}, if 1 G Q(£, A) then r(¥(0)%x,i) = 
£\,A,i> from which it follows that 

% Aj i G range(T) + VV, V£ G {1, . . . , fci} such that 1 G Q(£, A). 

If n\ t x > 1, let n be an integer in the range 1, . . . , n 1A — 1 such that for all m = 1, we 
have 

n tx ,m e range(T) + VV, V£ G {1, . . . , fci} such that m G Q(£, A). (6.13) 
We claim that this implies that 

%A, M +i e range(T) + VV, V£ G {1, . . . , fci} such that // + 1 G Q(£, A). (6.14) 
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1.1 



1,2 



"l,3 n l,4 



1.1 



1.2 



"l,3 D l,4 



1.1 



1,2 



1.3 



1.4 











* 

V*- 


• 

• * 

• '«Y 


• 


• 










• 

• • 


• 


• 


• 



+ 




Figure 4: Decomposition of Z G Mat cxc . In this example, k\ = 4. The elements which are 
not on the oblique segments are all zero. The first summand is T(Z). The second summand 
is in range(T). 

From this claim, and the fact that the above arguments are independent of the particular 
choice of A, it follows by induction that ( |6.12j ) holds. It thus remains to prove the claim. 

Suppose that £ G {1, . . . , ki} is such that /i + 1 G <2(£, A). Then a simple computation 
shows that 

T^(0)R^ +1 ) = n^+i + Gz >w , (6.15) 

where G^a^+i is a c x c matrix whose columns are all zero except possibly the (n^i + • • • + 
n i,x ~ £*) column, whose (x, u) component is given by the formula 



7x,«,C if m G V(x, A, /i + 1) 
otherwise, 



where the coefficients 7 XiU ,g are as in ( |6.10| ). Now, ^/(0)R^\ )fl+ i — T(^/(0)R^\^ + i) G range(T) 
implies that T(\l/(0)i?g i A, / i+i) G range(T) + W. From Remark 



x=i «e-p(x,A,M+i) 



(6.16) 
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from which it follows that ( |6.14| ) holds since E Xj \,u,fj,+i £ range(T) and fi x ,A, P ( A t+i,«) £ 
range(T) + W from the fact that p(fi + l,u) < ji + 1 and by the induction hypothesis. 



Proof of Proposition |6.2| The essential idea here is to decompose the proof into r 
separate blocks, where we use the arguments of the proof of Lemma ^T9] on each of the 
blocks. 

First, for all j G {1, . . . , r}, we let Wj denote the subspace of W consisting of matrices 
whose columns are all zero except the columns between iVj_i + 1 and Nj inclusively. 

As in the proof of Lemma |6.9| , we construct a basis { ^ij-^x, u ,m } of Wj, and we note that 
the space W in (|6.7| ) is equal to W\ © • • ■ © W r . From Remark |6/7|, for each j £ {1, . . . , r} 
we can choose vectors Vj^, . . . ,Vj^. which are such that the vector Uj(vj^) has a 1 in row £ 
and O's everywhere else. For each j G {1, . . . , r}, we then construct matrices Rj-^xm i n a 
similar manner as we did in Lemma |6.9| , whose only non-zero column is equal to one of the 
vectors Vj 7 \, . . . , Vj^ described above. We then define linear mappings Sj : Wj — ► ^(^(O)) 
as in ( |6.11|) , and it follows that Wj is isomorphic to Wj = Sj(Wj). 



Now, any Z £ Mat cxc can be partitioned as in 

Z 



i.e. 



Z 2 ,r 



(6.17) 



where Z VA is a (n P)1 + • • • + n p fcp ) x + • • • + n q> k q ) matrix. Rewrite Z as 



( z 1A o \ 

o z %2 



V 











( o z 1>2 z 1>r \ 

Z 2 ,l i?2,r 



+ 



/ 



y Z r \ Z, t 



r.2 







(6.18) 



/ 



where the second summand on the right-hand side of fl6.18| ) belongs to range(T) (Lemma 
6.4). Furthermore, following the same method as in Lemma for each j e {1, . . . ,r} we 



write Zjj as a sum of matrices as is illustrated in Figure |], where the second summand 
belongs to range(T). Thus, for Z as in (|6.17 ), we define T(Z) as the first summand on the 
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right-hand side of where the blocks Zj t j are of the form of the first summand in Figure 

f|. It follows that for any Z G Mat cxc , we have Z — T(Z) G range(T). 

Finally, we use an induction argument similar to that used in thejjroof of Lemma |T9] 
to show that for each j G {l,...,r}, we have Wj C range(T) + Wj. We thus define 
VV = VVi © • • • © VV r . This completes the proof. ■ 



7 Construction of a A-versal unfolding 



In this section, we will solve the problem which was posed at the beginning of Section || 
that is, given a linear homogeneous RFDE such as (|2.1|) and a set A of solutions to the 
characteristic equation (|2.2j) , find a A-versal unfolding ( |3.1| ) for (271). 



7.1 Main result 



The following useful result is proven in and 



Lemma 7.1 Let I be an interval in K and suppose that tpi, . . . , ip c G C(I, C ra ) are linearly 
independent functions. For each 9 a I, define <& c (9) = ((pi(9), ... ,(p c (9)). For a fixed tq G / , 
denote by q the rank of the n x c matrix $ c (tq). Then there exist c — q distinct points 
Ti, . . . , r c _q G J\ {t } swc/i i/iai t/ie n(c — q + 1) x c matrix co1($ c (to), . . . , $ c (r c _ g )) /ias ranfc 
c. 

Now, suppose (fx, . . . ,<p c are the basis elements of the space P in ( p.3|) , and let 3>(9) = 
((pi(9), . . . , <f c (9)), and I = [— r, 0]. Since $ = $5, it follows that the rows of $ are solutions 
to the vector ordinary differential equation £ = £ 5, and so the rank of the matrix $(6 I ) is 
independent of 6* JTJJ]. We now have the following 



Proposition 7.2 Let q = rank($(0)) ; and let r , . . . , r c __ g 6e as m Lemma \7.1\ . Then for 



any n x c matrix R, there exist matrices A. 



o- 



A,,..,, G Mat nxn such that 



c-q 



3=0 



Proof Define the linear mapping /C : Mat nX („( c _ g+ i)) — >■ Mat nxc by 



(7.1) 



/C(^l)=^-col($(r ),...,$(r c ^)). 
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If we associate Mat nx ( n ( c _ g+ i)) = C n ^ c q+1 ^ and Mat nxc = C nc in the standard way, then 
the nc x n(n(c — q + 1)) matrix representation of /C is given by 

/C~/ n «)(col($(To),...,$(T c _ g ))) T , 

whose rank is nc. Thus, JC is onto Mat nxc . For a given R £ Mat nxc , let .4 £ Mat„. X (n.(c-g+i)) 
be such that /C(.A) = R. The conclusion of the Proposition follows by partitioning the matrix 
A as A = (A , . . . , A c _ q ), where the Aj are n x n matrices, j — 0, . . . , c — q. I 



Remark 7.3 The crucial element in the proof of Proposition 7J: is the fact that the rank 
of the matrix col($(r ), . . . , $(r c _ g )) is equal to c. Lemma \7.1\ assures us that we can al- 
ways achieve this if we use q = rank($(0)) and if the delay times t q , . . . ,r c _ q are cho- 
sen appropriately. However, as is noted in j^j, in certain cases it is possible to achieve 
rank(col($(r ), . . . , $(r c _ 9 ))) = c with a value of q larger than rank($(0)). For example, 
it is possible to have the linear RFDE ($.J\) on C 2 such that P is ^-dimensional with basis 
matrix 



















where uj\ and u>2 are distinct real numbers. Thus, we have c 
However, for any r and T\ such that Tq ^T\, we have 



4 and rank($(0)) = 2. 



rank(col($(r ), $(ri))) = rank 





e" 




















gi<^2To 


e 


-IU12TQ 




e" 


-iuj\T\ 




















e 





4. 



From Proposition |6.2| , we know that 

Mat cxc = range(T) © VV, 

where W is isomorphic to W (see ( |6.7|) ), and the isomorphism 8 is described in the proof of 
Lemma |6i] and Proposition |6T2j Let . . . Q$} be a basis for W. From Proposition |7T2"| , 
there exist matrices A™, . . . , A™ £ Mat nxn such that 

c-q 

^(fi m ) = ^M/(0)A J m $(r J ), m=l,...,5, (7.2) 

i=o 
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where the set {£(Q m ) } spans a complement to range(T) in Mat cxc . Clearly the previous 
statement still holds if instead of being as previously defined, £ is any injective linear mapping 
from W into 71(^5/(0)) whose range is a complement to range(T) in Mat cxc . 

The following is our main result of this section, and follows from Theorem |5.1|, Proposi- 



tions |T2] and [7J| and the previous discussion. 

Theorem 7.4 Let {Q 1 , . . . , Q$} be a basis for W (see ( \6. Tj )), and let r , . . . , r c _ g be such 
that 

rank(col($(r ), . . . , $(r c _ g ))) = c. 
Let 8 be an injective linear mapping from W into ^(^(O)), such that 

Mat cxc = range(T) © £(W), 

(there exists at least one such £ ), and let A™ (m e {I, ... ,5}, j G {0, . . . , c — q}) be n x n 
matrices which solve t \7.£\) . For each m G {1, . . . , 5}, let L m be the bounded linear operator 
from C([— t, 0], C n ) into C n defined by 



c-q 



L m {z) = Y,Afz{r j ). 

j=0 



(7.3) 



Let C(a) be the S-parameter family of bounded linear operators from C([— r, 0], C n ) into C™ 
defined by 



C(a) = C + 2j a m L r , 



(7.4) 



m=l 



where the a m are complex parameters, and Co is as in Then j\3. i| j is a A-mini-versal 

unfolding of 



7.2 First order scalar equations 

In the case of first order scalar linear RFDEs, that is z t G C\ in ( [2.1| ), we are interested in 
the following question. Rewrite (|7.4|) as follows 

c-q / s \ 

C(a) = Co + E a ^ A T ( 7 - 5 ) 

j=0 \m=l / 

We wish to show that it is possible to find a change of coordinates which simplifies ( \l.5\) to 

c-q 
j=0 
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where f3j G C for all j = 0, . . . , c — q. Recall that B is a c x c matrix and 5 = dim(W). We 
have the following result. 

Proposition 7.5 // is a scalar equation, then 5 = c. 

Proof By Theorem 5.1 of ||, since n = 1, the number of Jordan blocks for each eigenvalue 
of B is 1. Thus c = n\ } \ + • • • + n r> \. Now, 5 is given by with kj = 1 for all j = 1, . . . , r. 
Hence the equality holds. I 



Theorem 7.6 Consider a A-versal unfolding of a scalar equation ( \2. 1\) given by ( \7.4) - Then 
there exists a change of coordinates in parameter space C s which brings (\7.5] ) to ( \7.Q ) where 
q = I, (3 = ((3 , . ■ ■ , Pc-x) with (3j G C for j = 0, . . . , c - 1 . 

Proof Since n — 1, q — rank($(0)) = 1 and therefore c delays To, ... , r c _i are necessary 
to solve equation ( |7.1|) . Let Ri, . . . , R$ be the 1 x c matrices such that S(fl m ) = ^/(0)R m . 
By Proposition |7.2| , 

c— 1 

R m = J2 A T^) 

3=0 

with AJ 1 G C. We can rewrite this equation as a matrix equation 

R T m = ($(r ) T , . . . , Hr c ^) T )A m = VA m . 

where A m = (A™, . . . , A^ 1 ) T and T is transposition. By choice of r , . . . , r c _i, the determi- 
nant of V is nonzero so that A m = V~ 1 R^ n . 

Let (3 = (j3 , . . . , Pc-i) T and a = («i, . . . , a c ) T . Since 5 = c, set 

P= (A\...,A c )a. 

Now (A 1 , . . . , A c ) = V^ 1 (Rj, • • • , Rj) is nonsingular since V~ x is a nonsingular c x c matrix 
and we claim that (Rj, . . . , -R^T) is also a nonsingular c x c matrix. Hence, this change of 
coordinates yields the result. 

We now prove the claim. Again by Theorem 5.1 of ||, the number of Jordan blocks for 
each eigenvalue of B is 1. The set of 1 x c vectors {Ri, . . . , Rs} is partitioned into subsets 
{Rj-,i,i,m\ Tn G <2(1, 1)}, for each j G {1, . . . ,r} and defined as in the proof of Lemma |Q| 
and Proposition |6.2j . That is, the only nonzero element in Rj-i^m is in the {n\^ + U2,i + 
. . . + rij t i — m + l) th column. Hence the vectors R±, . . . ,R$ are linearly independent since the 
unique nonzero element for each vector lies in a different column. ■ 
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7.3 Independence of unfolding on choice of basis 

Suppose we perform changes of bases in the spaces P and P* (see (|2.3|)), = U" 1 ^ and 
$# = $ f/ ; where U is an invertible c x c matrix. Then obviously we have (\l/*, $*) ra = 
and equation ( |4.10| ) transforms into 



x 



B*x + ^ # (0) [£(«) - Co) ($ # + h(a)U)x, (7.7) 



where 5* = C/ 1 BU . Let TJ/ be the linear invertible transformation of Mat cxc defined by 
Tu(M) = U- l MU, let T# : Mat cxc — ► Mat cxc be defined by T#(M) = [B#,M], and let 
71(^^(0)) be the set of all c x c matrices whose columns are in the range of the matrix 
^ # (0). Then it is easy to show that 

range(T # ) = TV (range (T)) 

and 

ft(^ # (0)) = Tu(K(y(0))). 
Consequently, from (|6.1|) , we have that 

Mat cxc = range (T # ) © VV # , 

where W # = TjjiVV). Let {ft±, . . . , Q$} be a basis for W (see (|6.7|)), and £ as in Theorem [FA 
Let A™, . . . , A™_ be n x n matrices such that 

c-q 

£(n m )= ^2^(0)A^{Tj), m = l,...,6. (7.8) 

3=0 

Then 



c—q c—q 



Tu(£(n m )) = U- l £(Sl m )U = U-^^Af^U = ^*(0)Af^*( rj ), m = 1, . . . , 6, 

3=0 j=0 

i.e. given the delay times To, ... , T c _ g which are such that Lemma |T. 1| holds, the matrices 
A™, . . . , A™ which solve ( |7.8| ) are such that the parametrized family ( |7.4| ) generates a A- 
mini-versal unfolding of (|2.1| ) independently of the choice of bases matrices $ and \& for P 
and P* respectively, provided that (\l/, $) n = I c . 
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7.4 Decomplexification 



In applications, it is usually the case that Cq in (|2.1|) and C{ot) in (|3.1|) are real, i.e. they 
are bounded linear operators from C([— r, 0], R n ) into R n . Although the previous theory has 
been carried out in complex spaces, it is straightforward to construct real versal unfoldings 
by a simple process of decomplexification of (|7.4|) . We need to assume that the set A defined 
in Section 2 is invariant under complex conjugation (one important example where this is 
always the case is center manifold reduction, in which case A is the set of all roots of (|2.2j ) 
with zero real parts). 

Theorem 7.7 Suppose that A = {A , A^, A^} where A is a subset of real eigenvalues and 
A h a subset of nonreal eigenvalues. Then, a real A-mini-versal unfolding of ([2. 1\) is given by 



where a p G R for p = 1, . . . , So, f3 S} f3 s+ s h G R for s = So + 1, . . . , So + d~h, L p is a bounded 
linear operator from C([— r, 0],R n ) into R n for p = 1,...,8q, and L s is a bounded linear 
operator from C([—t, 0], R n ) into C n , for s = 8q + 1, . . . , So + Sh- 

Proof The c x c matrix B can be decomposed as B = diag(-B°, B h , B h ) with c = co + 2ch, 
where B° is the cq x cq diagonal block of real eigenvalues Ao in real Jordan canonical form 
and B h is the Ch x Ch diagonal block of nonreal eigenvalues A^ in complex Jordan canonical 
form. 

We establish the following notation for the remainder of the proof. We let <pi, . . . , (p CQ and 
(pi, . . . ,ip Ch be bases for the generalized eigenspace corresponding to the eigenvalues of A 
and A h respectively, chosen so that the matrices $° = (<px, . . . , (p co ) and <& h = ((fi, . . . , ip Ch ) 
satisfy respectively $° = and $ ft = & h B h . Consequently, if we set $ = ($°, 

then the columns of $ form a basis for P, and we have $ = $ B. Now, let ip\ , . . . , ip* o be 
real linearly independent functions in P* (corresponding to Ao) and . . . be non- 
real linearly independent functions in P* (corresponding to Ah). If we denote \&°* = 
col(^, . . . ,ijj* o ) and = col(^*, . . . , V*J, then ^* = col(^ *, vj>^*) is a basis for P*. 
Define ^ = (\I/*, <I > )~ 1 \1/*, then (\l/, $) n = I c . Moreover, a simple computation shows that 
^ = col(\l/ , ty h , ty h ), where f is a Co x n real matrix corresponding to Ao whose rows are 
linearly independent, and ^ h is a Ch X n non-real matrix corresponding to A^ whose rows 
are linearly independent. 

We define three projections n , TEx, n 2 as in Remark |6w] where n corresponds to \1/ 
and has all real components while Hi and n 2 correspond to and ^ . Consequently, the 




(7.9) 



P =i 



s=6 +l 
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mappings 11! and n 2 are such that 

TT7 = n 2 . (7.10) 

Now let {fli, . . . ,Q$} be the elements of the basis of W used in the proof of Lemma 
|6.9| and Proposition |6.2| . Assume that these basis elements are ordered so that the set 
. . . , fl$ } corresponds to the block B°, {fls +ii ■ ■ ■ > ^s +s h } corresponds to the block B h , 
{Qs +i+s h , ■ ■ ■ , ^<5 +2<5 h } corresponds to the block B h , with 



Ch xc ®c h xc h ®c h xc h I P=1)---)O0j 

\ Cfe xc 0c h xc h 0c h xc h 



a. 



^ O co xco OcoXCft Ocgxc^ 

OchXco 0c h xc h ! ^s+6 h — "c h xco u c h xc h u c h xc h 

\ Ch xc c , jXCh Ch xc h / \ ChXCo Cfl Xc h O 



/ co xco C0 xc h CQ xc h 

Cfe xcn Oc.xc c ,x c , I , s = 6 + l, . . . ,5 + 5 h , 



where 0fc X £ is the k x £ zero matrix, is a Co x cq matrix with only one non-zero element, and 
is a Ch x 0^ matrix with only one non-zero element. Define R p = (R®, 0, 0) for p = 1, ... , 5 
where is the n x Ch zero matrix and R p is a n x Co matrix with only one non-zero column, 
corresponding to a vector t>z 6 M ra chosen as in Proposition |6.2| , such that \I/(0)-Rp has a 1 at 
the same position as the 1 in Q p . Note that vi can be chosen in ~$L n since n is real. Now, 



for all p — 1, . . . , 5 Q , by Proposition |7[2| we can find matrices A\, . . . , ^ such that 



c-q 



(ift 0, 0) = J] A$ ($ (r,), ^(r,), ^( Ti )), 

j=0 



where we have written $(tj) = ($°(r 7 ), $ ft, (r 7 ), ^^(rj)). The matrices A p - can be chosen to 
be real matrices for all p = 1, . . . , 5q. We then let L p be the bounded linear operator from 
C([-r, 0], W l ) into M n defined by L p (z) = £r« A v - z{r 3 ) for p = 1, . . . , 5 . 

Now, define i? s = (0, i?^, 0) and -R s +<5 h = (0, 0, R 1 ^) for s = 5 + 1, . . . , 5 + 5^, where the 
first is the n x Co zero matrix and the second designates the n x zero matrix. The 
nx Ch matrix R^ has only one non-zero column, corresponding to a vector v \ G C™ chosen as 
in Proposition |6.2| , such that \l/(0)i? s has a one at the same position as the 1 in Q s . It then 
follows (from (|7.10|) ) that ^(0)R S has a one at the same position as the 1 in £l a +s h - Then, we 
set £(fl s ) = ^/(0)R S and S(fl s+ s h ) = ^(O)^^. From Proposition [T2"| , there exist matrices 
Aq, . . . , A s c __ q such that 

(0,^,0) = g A^^(r J )^ h (r J ) 1 WH^)), 

j=0 
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which we rewrite as the system 

c—q c—q c—q 

E a j $0 ^) = °> E A i ^ = R *> E = o. 

3=0 j=0 j=0 

If we take complex conjugates of the above system, we get 

c—q c—q c—q 

E^°(^) = o, E^^)=^ E^ $/l ^) = > 

j=0 j=0 j=0 

which is equivalent to 

(0,0 ) ^) = E^($°(r,),$ ft (r J ),$^)). 

j=0 

Thus, 

c-q 

R s+s h = E AS ^)- 

3=0 

From this, we conclude that the matrices A s +5h , . . . , A s c ^ h of the decomposition 

R ^ h = E A ? 5h ^) 

3=0 

can be chosen so that A S j +5h = A?, s = 5q> + 1, . . . , So + d~h, j = 0, . . . , c — q. Now, for 
all s G {5o + 1, . . . ,So + 5h}, we let L s and L s+ s h be the bounded linear operators from 
C([-r,0],R n ) into C" defined by L s {z) = YT=o A i < T o ) and L s+Sh {z) = Yj=l z(r,). 
Taking real and imaginary parts of L s as in ( [7.9| ) yields a real A-mini-versal unfolding. ■ 

Corollary 7.8 If A = A 0; then a real A-versal unfolding of j\2.1\ ) is given by ( \7.Q) with 
f3 s = (3 S+Sh = for all s = 5 + 1, . . . ,S + 5 h . 

8 Examples 

In this section, we illustrate our theory with several examples. 
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Example 8.1 Consider the scalar delay differential equation 

x(t) = Co(x t ) = x(t) - x(t-l), (8.1] 



which has been studied extensively || [T^] . The characteristic equation has a double zero 
eigenvalue, therefore the center eigenspace P is two dimensional, i.e. c = 2. A basis for P is 
given by $ = (1,6*), and we notice that 

rank(col($(0), $(-1))) = rank ( * °\=2 = c 

By ( |6.6|) , 5 = dim(W) = 2, and if follows by Theorem |7.6| and Corollary |7.8| that a real 
A-mini-versal unfolding of (|8.1| ) is given by 

x{t) = £(a)(x t ) = x(t) — x(t — 1) + aix(t) + a 2 x{t — 1) 

where oti, a 2 € M. 



Example 8.2 Consider now the first order scalar equation 

x = C x t = A x x(t - n) + A 2 x(t - r 2 ) (8.2) 

where Ai, A 2 are real parameters and the delays n, r 2 are positive. Belair and Campbell 
have shown that (|8.2|) has points of nonresonant double Hopf bifurcation in parameter space 
(Ax, A 2 , Ti, r 2 ) with eigenvalues ±a>ii and ±a; 2 i. The center eigenspace P is four-dimensional, 
thus c = 4. 

Suppose that ( gg| ) has a double Hopf bifurcation at (A{, A%, t£ , r|), with rj* 7^ r 2 *. By 
equation ( |6.6|) we have that 5 = 4, and Theorem [7.6| implies that a complex A-mini-versal 
unfolding of ( |jj.2| ) has the form 

x = C(a)x t = (A\ + a>i)x(t — r*) + (A* 2 + a 2 )a;(t — r^) + a 3 x(t — r 3 ) + a 4 x(t — r 4 ) 

for suitable choice of r 3 and r 4 , where ctj G C for j = 1, 2, 3, 4. 



Example 8.3 The second order scalar delay equation 

x(t) + ax(t) + (3x(t) = f(x(t - r)), 



where / is a smooth function such that f(0) = 0, was studied in Longtin and Milton |T6 



as a model for the pupil light reflex. In Campbell and LeBlanc U, it was shown that the 
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linearization of this equation about the trivial equilibrium, which we write in first-order form 

as 

xi{t) = x 2 {t) 

%2(t) = —i3x 1 (t) — ax 2 (t)+Axi(t — T) 
has a double Hopf point with 1 : 2 resonance (eigenvalues ±i and ±2i) at parameter value 

5 -3 



(8.3) 



(a,/3,T,A) 



0, -, 71, 

' 2' ' 2 



•4) 



Using Theorem |5.1| , it is straightforward to show that if we treat a, ft, r and A as unfolding 
parameters which vary in a neighborhood of the point (|8.4j ), then ( |3.3|) generates a complex 
A-mini-versal unfolding for the singularity at parameter value ( |8.4| ). Proving this amounts 
to constructing the 20 x 16 matrix S in Theorem 5A and showing that its rank is 16. 

However, the goal of this example is to compute the real A-mini-versal unfolding for the 
singularity (|8.3| ) (at parameter value (|8.4Q ), which results from using the decomplexification 
procedure described in Section O. 



We have that a complex basis for the center subspace is given by 

JO „2i6 r ,-i0 „-2iQ 



e e 
ie ie 2ie 2ie 



e 
-it 



e 

-2te~ 2ie 



( 


e -is 


ie ls 


\ 




e -2is 


2ie~ 2is 






e is 


—ie ls 




\ 


e 2is 


-2ie 2is 


J 



A basis for the adjoint problem is given by 



which we renormalize by defining 

The result is such that 

/ -2i(3vr + 4i) (-3vr + 4i) 2 -32 n + 128 % + 6 tt 3 - 8 in 2 \ 

-i (-3 vr + 8 i) (3 7T + 8 if 64 tt - 16 m 2 + 64 % - 6 tt 3 

2i (3 tt - 4i) (-3tt -4i) 2 -32 tt - 128 % + 6 vr 3 + 8 m 2 



(0) = K 



\i(-3 



TT 



\i) (3tt - 8i) 5 



5.5) 



64 7r + 16i7T 2 -64i-67r 3 / 
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where k = (9tt a — 32tt 2 — 256) 1 . The matrix B is given by 



/ i 

2i 



\ o o 



\ 



-i 

-2i ) 



A basis for W is given by the following four matrices 



/ 1 








o\ 




( 








o \ 



















1 


















































U 








o) 




V o 








o J 


/o 








o\ 




( 








o \ 


















































1 



















\o 








oy 




V o 








1 J 



For any tq and t\ such that tq ^ Ti, the matrix col($(ro), $(ti)) has rank 4. Therefore we 



and T\ = — 7r. Following Remark pTFL we construct II,- : C 



choose, for example, r 
j = l,...,4: 

nj(v) = ^i(o)vi + tp j2 (o)v 2l j = 1, . . . , 4 

where ^(0) are the elements in the matrix \l/(0) in ( |S.5|) , and v = (t>i, t>2) T £ C 2 . 
we define 

1 




Therefore, 



Ri 



i-^n(o) 
^12(0) 



R2 



Rz 








1 



1-011(0) 



^12(0) 

A simple computation shows that 








R i 




Ri 




—1 

1-011 (0) 

012(0) 



$(0) + 



-1 

1-011 (0) 
Vl2(0) 



^22(0) 



1-011 (0) 
012(0) 



-It 
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and 



R 2 = o \ 2(l-^i(0)) o. l-fe(O) $ ( ) + o 2(1-^(0)) „. HMO) ) ^(-7T) 



Thus, if we set 



^22(0) ^22(0) / \ ^22(0) ' </>2 2 (0) 



l-^n(O) i-fe(o) 

Pi + I7i = ; — 777T — , P2 + «72 - 



iMo) ' " ' ^22(0) 

then the operators L m in ( |T.3|) are 

Li(z) = i [ 1 ~ l | z(0) + ~ | _1 I :(-tt). 

4^1 + ^71 7i-^i 7 4\-Pi-iji - ll +if3 l 1 

L 2 (z) = -( 2 ~ % | z(0) + i I 'i :(-tt). 

8 ^ 2/? 2 + 2z 72 2 72 - 2i(3 2 ) %\2(3 2 + 2i l2 2 l2 - 2if3 2 ' 

Ls — Li, L 4 = L 2 . 

Thus, the decomplexification procedure yields the following real A-mini-versal unfolding of 
the singularity ( |B.3|) (at parameter value ( |S.4|) ) 



+ 



xt(t) \ f x 2 {t) 

X 2 {t) J \-\x 1 {t)-\xx{t-'K) 



1 o o -i UW-.H 
A 71 y V ^ JJ\ ^(t) - x 2 (t - n) 1 

2 \ / -1 V, xW+xiCt-.) 
2(3 2 2 72 / 1 2 72 -2(3 2 jj \ x 2 (t) + x 2 (t - n) 



Appendix 



A Proof of Lemma |6t4 



We use the following inner product on the space Mat CX c- 

(Mi, M 2 ) = trace(M 1 M*). (A.l) 
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A simple computation shows that with respect to ( |A.1| ), the conjugate transpose of T is 
given by T*(M) = [B*,M]. From the Fredholm alternative, we have range(T) = ker(T*) ± . 
Recall that ker(T*) has been characterized in Lemma |6l| So Y G ker(T*) ± , if and only if 
Y is orthogonal to all elements in a basis of ker(T*). We construct a basis for ker(T*) as 
follows: for j G {1, . . . , r}, denote 



n (j) = number of distinct oblique segments in the block Aij, as in Figure 



and give some ordering to these oblique segments, numbering them from 1 to n Q (j). Then, 
for j G {1, . . . , r} and £ G {1, . . . , n Q (j)}, we define to be the matrix structured as in 
( p.3| ), and such that all diagonal blocks are zero except the block j; and in this block j, the 
only oblique segment which is non-zero is the £ th segment, whose elements all have value 1. 
Then {Mj/} forms a basis for ker(T*). 

Now consider Y G Mat cxc , which we partition as in (O). Multiply Y on the right by 



M* f} for some j G {1, . . . , r}, £ G {1, . 
is such that 

/ 




YM te 



n (j)}. Then YM* e is partitioned as in ( |6.4| ), and 












V o 



• ° yrjM* j)t 

T ' 

jth vertical block 
in partitioning 



0/ 



where Aij/ is the j diagonal block in the decomposition (|6] 
trace(FM^) = trace (yj,jM* u ). It follows that Y G ker(T* 



of M jt £. Thus, (Y,M j>e ) = 
- if and only if for all j G 



{!,..., r}, the block in the partition 



is such that trace (yj,jAi%) = for all 
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G {1, . . . , n (j)}. Using the partition illustrated in Figure |I|, we write 



y 



^3,3 
^3,3 



3? 



1,2 

3 
2.2 



-yl.fcj 

^3,3 

-y2,ftj 

^3,3 



\ 



\yj;. 



^3,3 



(A.2) 



and 



A4 



1,2 

2,2 



■ykj,kj 



M k /f j 



(A.3) 



/ 



where 3^f ',■ and A4 



<P>9 



are Tij^p x Tij qi 



for p,q & {1, . . . , fcj}. Let p and g be such that AWj'| is the 
unique block in (|A.3| ) which has the non-zero oblique segment. Then yj,jM.*j £ is partitioned 
as in ( |A.2[ ) and has the form 



/o 










y)j ■ (M* 



y 



2,9 



(^0 












V • ■ ■ y jf ■ C^??)* 

T 

pth veT iical block in 
partitioning ( |A.2 ) 



It follows that Y G ker(T*) _L if and only if Y is of the form 



and for all j G {1, . . . , r} 



and for all £ G {1, . . . , n (j)}, we have trace(3^^ • (A4^'|)*) = 0, where p and q are such that 
A4y'| is the unique block in ( |A.3j ) which has the non-zero oblique segment. Now, suppose the 
block M. P a'1 is such that the oblique segment whose elements are all l's is the one which passes 
through the v th leftmost element in the bottom row of the block, where v G {1, . . . , n,j t q}. 
Then a simple computation shows that y™ ■ (•M^'J)* is an rij^ x rij^ matrix whose first 



36 



rij t p — v rows are zero, and the remaining rows are the rows 1, . . . , v of yf] ■ It follows that 
trace(>'Jj ? • (M^j)*) is the sum of the elements in the oblique segment of the block y?'J 
which is the same oblique segment as the non-zero segment of M-^'j, from which we get the 
conclusion. I 
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